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Abstract 

In most proposals for the generation of entanglement in large ensembles of atoms via projective 



Cn ■ measurements, the interaction with the vacuum is responsible for both the generation of the signal 



that is detected and the spin depolarization or decoherence. In consequence, one has to usually work 
in a regime where the information aquisition via detection is sufficiently slow (weak measurement 



,. . regime) such as not to strongly disturb the system. We propose here a four-wave mixing scheme 

\^ . where, owing to the pumping of the atomic system into a dark state, the polarization of the 



ensemble is not critically affected by spontaneous emission, thus allowing one to work in a strong 
measurement regime. 



I. INTRODUCTION 

Considerable attention has been given in recent years to the generation of multi-particle 
entanglement in large ensembles of long-lived atomic spins. Some authors have proposed 
to achieve this goal using interactions between atoms and light, where a quantum state 
excha^ge can ta.e place fiflSflQ. .ead,„g to the preparation of a des.ed eoUect.ve spin 

state. A different set of proposals make use of an appropriate measurement on one of the 



field's observables that leads to the collapse of the ensemble onto the desired 



ent ang 



Ifl 



11 



ed spin 



13 and 



state. In general, these schemes fall into two categories: conditional 

deternnnisfc fl U Q. In a conditional entanglement generation schcne. the prepared 

state of the atoms is conditioned by the outcome of the measurement on the field state. The 
random character of the state resulting from the measurement back-action can be removed 
if one performs a continuous quantum nondemolition measurement of a spin observable and 
adding a feedback loop for a posteriori quantum state correction based on the detection 
outcome. In this way a deterministic (unconditional) quantum state preparation scheme 
can be realized, where the uncertainty in the final state is removed. 

In general, multilevel atoms with long-lived ground substates are used, where two of the 
ground sublevels form an effective two-level atom. The manipulation of the collective atomic 
spin, obtained by summation over the individual spins associated with each atom in the 
medium, is achieved by driving ground-excited state transitions using classical or quantized 
radiation fields. As a result of the interactions, a signal field is generated that reflects some 
quantum mechanical fluctuations in the atomic ensemble; its detection can give information 
about the atomic ensemble state. In the off-resonance regime, an effective Hamiltonian can 
be found with a coupling between atoms and signal field that is proportional to the strength 
of the atom's coupling to the vacuum and also to the amplitude of the driving field. An 
increase in the control field's amplitude, therefore, would seem to allow one to generate 
optimal entanglement. However, the downside of using coupling through excited levels is 
that spontaneous emission comes into play, leading to a rapid decoherence of the system. 
In consequence, it is necessary to limit spontaneous emission losses to a small value, which 
forces one to work in a regime of weak coupling of the atomic system to the field system to 
be measured, and only weak entanglement can be obtained. 

The competition between spontaneous emission and measurement strength is best illus- 



trated in the case of spin squeezed state [16[ generation. The challenge there is to reduce 
fluctuations in a spin component orthogonal to the mean spin, while keeping the average 
spin length large. The measurement strength controls the reduction of fluctuations, while 
spontaneous emission leads to a diminishing of the spin length. In a recent publication |l7|. 
it has been shown that, for the case of a pencil-shaped medium with Fresnel number close 
to unity, in the regime of small decay, optimal results are limited by the resonant optical 
depth of the sample. 

Despite the interplay between measurement strength and spontaneous decay, it is pos- 
sible to imagine a situation in which spontaneous decay does not limit the value of the 
measurement coupling strength. In one such scenario, the system is prepared in a dark state 
which is preserved during the interaction by a convenient choice of driving fields. This is 
the situation presented in this paper, where a quantized signal field, generated via four- 
wave mixing in a double A atomic system in a pencil-shaped medium, is entangled with the 
collective atomic state. The generated signal pulse reflects fluctuations in the population 
difference between the two ground substates. The measurement on the signal field photon 
number can, therefore, give information on the z component of the atomic spin, projecting 
the system into either a spin squeezed or Schrodinger cat state. At the same time, as op- 



posed to the situation illustrated in 



3, 



the collective x polarized spin state in which the 



system is initially prepared, is very nearly a dark state for this combination of fields and 
is subject to minimal decay. No severe limitations on the measurement coupling strength 
are therefore necessary, and results similar to the ideal case presented in |l2| are obtained. 
With the assumption of perfect detection, the Heisenberg limit is the ultimate limitation to 
the squeezing parameter. 

The paper is organized as follows: in Sec. II the proposed scheme is described and 
an analytical expression for the atomic operator giving rise to the signal field is obtained; 
spontaneous emission effects are also discussed. In Sec. Ill, an expression for the signal field 
amplitude operator is obtained, which is shown to reflect atomic population fluctuations. 
An effective Hamiltonian necessary for a wave function description of the problem is derived 
in Sec. IV, while in Sec. V the generation of entanglement via conditional measurement of 
the signal field is discussed. Some conclusions are presented in Sec. VI. 



II. SCHEME AND METHOD 

A pencil- shaped atomic medium aligned along the z axis (left end situated at 2; = 0), with 
transverse area A, length L, and density n^ is considered. The internal structure of an atom 
[Fig. n^b)] is a double A scheme with ground levels 1 and 2 and excited levels 3 and 4. Three 
classical laser pulses having duration T ^ L/c (where c is the speed of light), wave vectors 
ki, k2 (pumps) and kp (probe), and frequencies fii, ^2 and fip are simultaneously shined on 
the atoms. We consider an off-resonant regime with one-photon detuning A for the pump 
fields, one-photon detuning Ap for the probe field, and two-photon Raman detuning 5 [see 
Fig. n^b)]. The three classical waves mix inside the atomic medium to generate a multitude 
of secondary waves; the one which propagates along the positive z direction [see Fig. lU^a)] 
and which can be viewed as a reflection of the first pump wave off a spatial grating produced 
by the second pump and the probe is of interest and denoted as the signal wave. This wave 
is radiated on the 2—^3 transition, and has frequency fi^ = i7p — (fii — ^22) and phase 
matched wave vector k^ = (i^^/c) z = kp — (ki — k2). A few methods can be employed to 
separate the signal wave from the probe and pumps. The simplest one [illustrated in Fig. 
^a)] requires an increase of the angles made by the three primary wave with the z axis 
sufficient to provide a clear angular resolution. The other two involve both the use of a 
polarization beam splitter (to distinguish between probe, pump 1 and signal) or a spectral 
filter (that can distinguish between pump 2 and signal). Finally, a photodetector (PD) is 
used to detect the photon number of the signal field. In the following, we consider the 
dynamics of an atom located inside the sample, at position r. 

The three incoming pulses are represented as c number traveling waves 

E«(r,t) = - [Ei(r,t)e^(^'-'-^»*) + E*(r, t)e-^(''»-"'^"*)] q, (1) 

for i = 1, 2,p, with polarization unit vectors q. The matrix elements of the atomic dipole d 
between atomic states j and j' (along q) is denoted by d^y = {j |d-'ej| j'). The interaction 
between fields and atom is described by the (spatially dependent) Rabi frequencies 

Ai(r,t) = ^^k^^e^i^'- = xi(r,t)e^'^»-, (2a) 

A2(r, t) = ^k^^e^k.r ^ ^^(^^ ,)^.k..r^ (2b) 
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FIG. 1: (a) The three classical pulses indexed by 1, 2 and p are incident on the sample at an angle 
giving rise to a signal field that propagates along the cylinder's axis. Detection of photon number 
takes place at the photodetector (PD). (b) Illustration of the internal structure of an atom as a 
double A scheme, (c) The signal field is respresented as a reflection of the probe [p) off a spatial 
grating generated by the pumps (1 and 2) 



Assuming that the pump and probe pulses are long [cT ^ L] and wide (transverse area 
larger than A), we can neglect the slow variation with r in Xj(r,t); the field envelope is 
replaced with its value at r = 0, and Xj(0,t) is replaced with Xj{t)- 

Atomic operators are defined as an = \i) {i\ (population operators), for i = 1,2, 3, 4, and 
'^ij = K) 01 (coherences), for i,j = 1,2,3,4 but with i ^ j. Ignoring the coupling to the 
vacuum for the moment, the Hamiltonian that describes the evolution of the atom driven by 
the three fields is a sum of the free Hamiltonian {Hq) and the classical field-atom interaction 
Hamiltonian (Vc) 

Hq = fUjJ2l<y22 + ^310^33 + ^4lO"44, (3a) 

K = -n [Ai(r, t)e-*^^V4i + At(r, t)e'^^Vi4] - (3b) 

- h [A2(r, t)e-'^^V42 + A;(r, t)e''''''a24\ - 

- h [Ap(r, t)e-^^-V3i + A;(r, t)e^^-Vi3] • 

In the Heisenberg picture, the rapid time variation in the coherences is removed: (T14 = 



0-146 *^l*, 0-24 



0-246 '^^\ 0-13 



0-136"*^?'*, 0-12 = 0-i26-*(^i-^2)t^ cr34 = cr346-*(^i-^f')* and 



Cr23 = 0"236^ 



-i[Q.p-{ni-n2)\t 



Equations of motion for the slowly varying atomic operators (o-jj) 



are obtained (dropping the tildes) as 
d 



d 

d 

d 

d , 

d 

d 

d . 

d , 

d 



-zAi(r, t)o-4i + zAt(r, t)o-i4 - zAp(r, t)o-3i + iA*(r, t)o-i3, 

-■iA2(r, t)o-42 + ^^*2{r, t)o-24, 

iAp(r,t)a3i - zA*(r,t)o-i3, 

iAi(r,t)o-4i -iA^(r, t)o-i4 + zA2(r, t)o-42 - ■iA2(r, t)o-24, 

-zAcri4 +iAi(r,t) [o-ii -0-44] +zA2(r,t)o-i2 - zAp(r, t)o-34, 

-i (A + 5) CT24 + iA2(r, t) [0-22 - 0-44] + iAi(r, t)o-2i, 

-zApO-13 + zAp(r,t) [0-11 - 0-33] - zAi(r,t)o-43, 

-i (A - Ap) 0-34 + iAi(r, t)o-3i + iA2(r, t)o-32 - iA*(r, t)o-i4, 

-i (Ap + 5) 0-23 + «Ap(r, t)o-2i - iA2(r, t)or43, 

i5ai2 - iAi(r, t)a42 + iA;(r, t)ai4 - 2Ap(r, t)a32. 



(4a) 
(4b) 
(4c) 
(4d) 
(4e) 
(4f) 

(4g) 
(4h) 

(4i) 
(4j) 



We are interested in the time and spatial modulation of the a23 coherence, which is re- 
sponsible with the generation of the signal. The rapid time variation of 0-23, at frequency 
VLs = ^p — (fii — ^2) has already been eliminated. The terms giving a spatial modulation 
at the correct, phase-matched wave vector k^ = kp — (ki — k2), have to be identified in Eq. 

(Hi))- 

We proceed to describe the dynamics of the atom by using a perturbative approach to 

solve Eqs. (JH), in which terms up to the third order in xl ^i ix = Xi5X2,X3 and Aj = 

A, Ap, 5) are kept. We start with an x polarized state of the two-level system formed by 

the ground sublevels 1 and 2, where o-ii(O) = aii = 1/2, 0-22(0) = a22 = 1/2, 0-12(0) = 

0-^2 = 0-21(0) = o-gi = —1/2. Owing to the assumptions of large detunings and long duration 

pulses, the upper state population operators an can be neglected. The coherences between 

ground and excited states adiabatically follow the fields. Assuming 5 ^ A, Ap, the first 



order solution of Eqs. (jl^,f,g,i) are given by: 



aSlfi . ^.!. + ^.L (5a) 

4'W==^.5, + ^.,°. (5b) 

,(.)(,) ^ ^„l, (5c) 

4l\t) - ^<. (5d) 

To second order, as observed before, the excited state populations have derivatives which 
are identically zero; however, the ground states rate equations are 






At(r,t)A2(r,t) ^ Ai(r, t)A^(r, t) „ 

■CTio : 0"n 



12 A 21 



A '" A 



(6a) 



|4?W^-|.!?W. (6b) 

An important observation can be made at this point. The system formed by the two 
ground sub levels is driven by an effective field with a Rabi frequency [xl(t)x2it)/^] niulti- 
plied by a spatial phase, dependent on the atom's position inside the medium: e-«{ki-k2)r_ 
In the copropagating pumps geometry, this spatial modulation is negligibly small over the 
length of the medium since |ki — k2|L ~ uj2iL/c ^ 1. In addition, assuming Xiif) and 

• (2) 

X2{t) are real, the resulting effective Rabi frequency is real. The a^ is, in consequence, 
vanishing for any atom inside the medium, independent on its location. Also, owing to the 
fact that equal the ground state populations are equal, the cri2vanishes to this order as well; 
in consequence the system stays in a dark state with equal populations and coherence along 
the X axis. 

The two coherences, (734 and (T12, that act as sources in Eq. (jH), are not driven directly 
by the field; their change from initial values is a second order contribution. They can be 
derived from Eqs. (jlhj), which give 



<y?ht) ^ -4l^ [Ai(r,t)a?i + A2(r,t)a?2] , (7a) 



'■p 



J^h+\ ,0 , Ai(r,t)A*(r,t) ^ p .. 

^12 {t> = ^12 + -^ [^11 - ^22] • ( 'bj 

We are now in position to evaluate the third order approximation of cr23- Replacing Eqs. 



(jTj) in Eq. (jU), one obtains 

(3) Ap(r,t) At(r,t)A2(r,t)Ap(r,t) . ^ „ . 
c^23 (^) ^ ^^^21 + ^^^^ kn - ^22j - (8) 

Al{r,t)A2{r,t)Ap{r,t) q [A^lr, t)|' Ap(r, t) „ 

Two of the terms in the above expression [first and fourth in the right hand side of Eq. (jHJ] 
describe a field propagating in the direction of kp, which is not phase matched [kp ^ fis/c]. 
Both other terms give rise to a phase matched signal field; however, in the limit b ^ Ap 
the third term is negligible compared to the second one, and is dropped. With the notation 
Oz = (o'22 — o"ii) /2, the expression of (T23 can be simplified 

At(r,t)A2(r,t)Ap(r,t) ^_^ 

^23^ = -^^ [fTii-a22j- (9) 

The analysis is not complete before the role of spontaneous decay is properly identified. 
The conditions imposed on the fields guarantee that the effective driving pulse doesn't 
remove the system from the dark state; spontaneous decay can still destroy the coherence, 
as is the case in [l7|. However, we start by making an observation on a simple system 
of a A atom driven by two equal fields, in which, on each of the two transitions, Raman 
and Rayleigh scatterings cancel each other, leading to a state not affected by spontaneous 
emission. Our case is similar to this, although not completely identical. In the absence 
of the probe field, an initially balanced state (equal populations in the ground substates) 
with —1/2 coherence, would be preserved by choosing equal amplitude pumps. The probe 
field provides an imbalance in the system, which can be compensated by choosing a field 
strength on the 2^4 transition larger than the one on the 1 — > 4 transition by an amount 
that cancels the effect of the probe. This condition is 7'(xi — xi)/^^ = 7Xp(^)/^p ^^^ is 
obtained by imposing the condition that the decay terms in the rate equations for crn{t) and 
(722 (i) are identically zero. In addition, the coherence between the ground sublevels follows 
the fields with a slowly varying value of —{xi/x2)/2. By limiting the intensity of the probe 
field to small values compared to the pump field intensities, this coherence stays close to the 
maximal value of —1/2 throughout the interaction. The maintainance of coherence is the key 
feature of this level scheme, allowing for much better spin squeezing than in other projection 
schemes. A similar idea was used to improve spin squeezing in cavity-field interactions [l8| 



III. EMITTED FIELD 



A wave equation for the signal field can be written where the polarization of the medium 
acts as a source. Defining the positive frequency part of the polarization resulting from a 
single atom (denoted by a) located at position r as Pa (t) = Pa{t)e^^^^^~^^''\ its envelope 
is given by 

Pa{t) = h[d2^a2z{t) + h.c] (10) 

The assumption of Fresnel number close to unity for the pencil-shaped medium can be 
invoked now; this leads to a one- dimensional behavior of the propagation of the signal field. 
In consequence, as in Refs. |l9|, |20|, z dependent continuous operators can be defined by 
performing an average over infinitesimal slices in the transverse direction of the medium. 
The continuous polarization operator is thus defined as 

1 



P{z,t) = lim 



E Pa{t), 



(11) 



where AVz = AAz is the volume of a slice and the sum is performed over all atoms in the 
slice (number of atoms in a slice N^ = UaAAz). Continuous atomic operators can also be 
defined as 

0(^,t)= hm-L ^ dait). (12) 

Replacing the expression for cr23 previously derived in Eq. (jTUI) and making use of continuous 
atomic operators, Eq. (fTTj) becomes: 



P{z,t) = huadisfit) [au{z) - (y22iz)] . 



(13) 



where the notation f{t) = Xi{t)X2{t)Xp{t) I AAp6 has been made. As in Ref. [jjl, |2i|, the 
quantized signal field amplitude (the positive frequency part) can be written as E^^\z,t) = 
SsEs{z,t)e^^''''^~^''^\ where Sg = \/firtj2t^AL and Es{z,t) is a slowly varying envelope 
operator. The wave equation in terms of slowly varying field and polarization envelope 
operators can be written as 



d, + -dt 
c 



Es{z,t) 



2eo 



P{z,t). 



(14) 



Setting K(t) = hksd23f{t)/2eQ, we find that population fluctuations at each point in the 
medium are connected to the signal field by 



d, + -dt 



Es{z,t) = -iUaKit) [0-11(2;) - 0-22(2:)] . 



(15) 



If collective atomic operators are defined 

L 

Sz = ^j dz[aii{z)-(j22{z)], (16) 



the solution for Eq. (fT3j) is found [for derivation see Appendix A] to be 

E,(L,t)=E,(0,t)-z^^5,. (17) 

The result states that the signal field amplitude exiting the sample is amplified by a 
quantity proportional to the collective population operator. The factor 2[K{t)/A\Sz can 
also be reexpressed as na^Sz/Na\L, which shows a linear increase with the length and 
atomic number density of the sample. The intensity of the emitted field can be calculated 
as an expectation value 

%{m)El{L,t)) = ^-^^^{Sl). (18) 



An initial population imbalance in the ground substates gives rise to a signal quadratic 
in the number of atoms in the sample. However, when the collective state of the system is a 
coherent one, polarized along the x direction, for example, population fluctuations only are 
reflected in the emitted field. The variance of Sz is in this case iVf /4, which leads to a gain 
in the field intensity, linear in number of atoms. 

IV. EFFECTIVE HAMILTONIAN AND ENTANGLEMENT GENERATION 

The conditional atomic generation process is similar to an EPR-type experiment, where 
entanglement is created between two subsystems (medium and signal field, in our case) 
by means of an interaction that lasts a finite time (duration of pulses); a measurement 
(detection of signal photon number) is performed on one of the subsystems (field) long after 
the interaction has ceased. Consequently, the other subsystem (atoms) is projected onto 
the state entangled with the state indicated by the detection outcome. A wave function 
approach (or density matrix, when imperfect detection is accounted for) can be taken to 
describe the coupled evolution of the two subsystems, while a continuous measurement 
theory is particularized to this case to describe the nondeterministic evolution of the system 
during the detection process. 

10 



The results of Ref. [17| are used in what follows. The complete details of the derivation 



of an effective Hamiltonian are found in Appendix C of Ref. 



17| , where a similar calculation 



is described. The measurement process, both under the assumption of perfect detection and 



including imperfect detectors, is also presented in details in Sec. V. of Ref. [17[. We are 
concerned here, rather with the main differences between our proposed scheme and the ones 
proposed elsewhere. 

The interaction between the signal field amplitude operator and the atoms in the sample 
is written in the Heisenberg picture. An integration over the transverse wave vector com- 
ponents of the generated field (allowed by the assumption of Fresnel number close to unity) 
followed by one over x and y leads to a one-dimensional formulation of the problem, where 
the continuous atomic operators are specified only by their z spatial location, while the field 
. las a transverse spatial extent A (matching the cross-sectional area of the medium) [see also 
20|]. The coherence between levels 2 and 3 is thereafter replaced, using Eq. Q, to lead to 
an effective Hamiltonian 

H.ffit) = hhit) (/dA:,4(A;,)e^(-^-^)*) S, + /i.c, (19) 

where h{t) = 4:Trd23£s/hy/A f{t), while dUk^) is a one-dimensional field operator defined 

1 



as 



d\{h 



--Jdxdyjdk^dkyax{k)e '"''e 



tUxX IfZyXf 



(20) 



27rv^A 
and satisfying the following commutation relations d\{kz)., d\,{k'^) = 5{kz — k'^)6x\'- 

With the observation that the Hamiltonian commutes with itself at different times 
[Heff(t),Heff{t')] = 0, the evolution operator over the duration of the interaction can be 
expressed in a simple form: 

o T 



U{T) = exp 



-TjdtH^ffit) 



(21) 



The time integral brings the Fourier components of f{t) the incident pulse field envelope 

T 

J dte^^^''~^^^ f (0 , t) ~ F{ijJk — f^). The integral over k^ in Eq. (J19p can be now represented by 



an effective one-photon creation operator with carrier frequency Vt^ and duration (cA/c)^ ~ 
T defined as 



.1/2 



Wfit)\ 



-.jdk,F{ujk-n)dl{k,), 



(22) 



11 



and obeying the usual commutation relation [cy, cj] = 1. This leads to a simple form for the 
evolution operator 

U{T) = eM-^C{cl-Cy)S,]. (23) 

with 

1/2 

C = [47rd23^./?iVIc^/'] Cjdt \f{t)f\ . (24) 

The atoms-signal field system starts in an initial state with state vector 

s 

M=-S 

where the index a denotes states of the atoms, while the index / denotes states of the field. 
The initial state of the atoms is an eigenstate of Sx [operator which is defined similarly to Eq. 
^] with binomial coefficients A{S, M) = ^^{2S)\/{S + M)\{S ~ M)\ (where S = Na/2). 
After a period of coherent evolution governed by the evolution operator U(T) [Eq. pHjl ]. a 
collapse induced by a measurement with an outcome of n^ photons leads to the following 
state vector for the atoms (assuming 100% detection efficiency): 

WVJ. ^ E , ^(^■^OC^-'^)"--'-"'^- I I..M). (25) 

J2 |A(5,M)|'(CM)2"-e-(^^^)' 

'=-S 

This state vector describes spin squeezed states for rim = and Schrodinger cat states 




for Um > 0. The results are similar to the ones presented in 



til 



(the reader is refered to 



that publication for relevant discussions and graphs), with one important exception. In 



0, 



the value of C is limited to small values {C < ^y{naX'^L)/2Na) to insure that spontaneous 
emission does not lead to mean spin depolarization. Here, that restriction does not apply 
since the total coherence is not seriously affected by spontaneous emission. In consequence, 
the results obtained here overlap with the ideal case presented in [17'], where spin squeezing 
close to the Heisenberg limit and well-resolved Schrodinger cat states can be obtained. 

V. CONCLUSIONS 

We presented a probabilistic scheme in which the detection of photon number induces 
the collapse of the quantum state of a collection of atoms onto either a spin squeezed or a 
Schrodinger cat state. The main result of the paper is that spontaneous decay does not play 

12 



a critical role in the decoherence of the system, therefore allowing one to obtain a squeezing 
parameter close to the Heisenberg limit (assuming perfect detection). This has been done by 
maintaining the system in a dark state by choosing the appropriate configuration of driving 
pulses. 
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VII. APPENDIX : SOLUTION FOR THE FIELD EQUATION 

Using the Laplace transform with respect with z (defined as o(s) = £{0(2;)} = 
JdzO{z)e^'^'') and defining e{s,t) = J^ \ Es{z,t) > and ct^(s) = C {az{z)} Eq. (fT3j) becomes 



(Al) 



dtCs^s, t) + cses{s, t) = cEs{0, t) — icK{t)az{s). 



A formal integration results in 



e,(s, t) = e,(s, 0)e-"'^* + / dt' e-""'^'-''^ cEsiO, t') - icK{t')a,{s) 





(A2) 



The time dependent term K{t') contains the slow varying field envelope, which is evaluated 
at time t leading to 



1 



e,(s, t) = e,(s, 0)e-'=^* + - [l - e-'^'] E,(0, t) - tK{t)a,{s 



(A3) 



Applying the inverse Laplace transform £ ^, a general solution of Eq. (|15|). for arbitrary z 
and t, is found 



Ei+\z,t) = Ei+\z - ct,t)h{z - ct) + ^i+^(0, t) [h{z) - h{z - ct)] 



(A4) 



iK{t) 



z—ct 



dz'ajz') - < / dz'a^iz') } h{z - ct) 

K Q 

We evaluate now the field at the sample exit {z = L). In the limit of long pulses (T ^ L/c), 
the Heaviside function h{L — ct) is zero for most of the time and will therefore be ignored. 



13 



In terms of the collective operator defined in Eq. ()16p . the field becomes: 

Ei^\L,t)=Ei^\0,t)-^^-^S.. (A5) 
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